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00 , Abstract: We present a construction of a 2)-dimensional Ricci-flat metric correspond- 

. ing to a (c? + l)-dimensional relativistic fluid, representing holographically the hydrodynamic 

regime of a (putative) dual theory. We show how to obtain the metric to arbitrarily high order 
using a relativistic gradient expansion, and explicitly carry out the computation to second 
I order. The fluid has zero energy density in equilibrium, which implies incompressibility at 

flrst order in gradients, and its stress tensor (both at and away from equilibrium) satisfies 
a quadratic constraint, which determines its energy density away from equilibrium. The en- 
^ ■ tire dynamics to second order is encoded in one first order and six second order transport 

^ . coefficients, which we compute. We classify entropy currents with non-negative divergence at 

second order in relativistic gradients. We then verify that the entropy current obtained by 
pulling back to the fluid surface the area form at the null horizon indeed has a non-negative 
divergence. We show that there are distinct near-horizon scaling limits that are equivalent 
either to the relativistic gradient expansion we discuss here, or to the non-relativistic ex- 
pansion associated with the Navier-Stokes equations discussed in previous works. The latter 
expansion may be recovered from the present relativistic expansion upon taking a specific 
non-relativistic limit. 
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1 Introduction and summary of results 

During the last year an interesting "holographic" connection between Ricci-flat metrics and 
fluids was uncovered: in [1] an approximate regular {d + 2)-dimensional Ricci-flat metric 
corresponding to a solution of the incompressible non-relativistic Navier-Stokes equations 
in (d -|- 1) dimensions was presented, valid to leading non-trivial order in a non-relativistic 
hydrodynamic expansion, and in [2] we provided a systematic and unique construction of 
this metric to all orders. The construction has been extended to first non-trivial order to 
spherical horizons in vacuum gravity [3-5], to de Sitter horizons [6] and to higher-derivative 
theories coupled to matter [7, 8]. Important earlier related works include the membrane 
paradigm [9-11] and the more recent construction of solutions of AdS gravity describing the 
hydrodynamic regime of CFTs in the context of the AdS/CFT correspondence [12-15], see also 
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Figure 1. The analogue of the thermal state in 
our construction is Rindler spacetime, with past 
and future horizons and "H^ respectively. The 
dual fluid lives on a constant acceleration surface 
Ec with flat induced metric. Lines of constant r 
and constant r in the coordinate system (2.1) are 
shown in grey. 



earlier results in [16]. Another instance of the fluid/gravity correspondence in asymptotically 
flat spacetimes can be found in the blackfold approach [17, 18]. Further developments were 
reported in [19-28]. 

The existence of fluid solutions in gravity is expected/predicted by holography on general 
grounds. A generic feature of QFTs is the existence of a hydrodynamic description capturing 
the long-wavelength behaviour near to thermal equilibrium. One then expects to find the same 
feature on the dual gravitational side, i.e., there should exist a bulk solution corresponding to 
the thermal state, and nearby solutions corresponding to the hydrodynamic regime. Global 
solutions corresponding to near-equilibrium configurations should be well approximated by 
the solutions describing the hydrodynamic regime at sufficiently long distances and late times. 
This picture in indeed beautifully realised in the AdS/CFT correspondence, where the thermal 
state corresponds to a bulk black hole [29] , and nearby solutions describing the hydrodynamic 
regime, corresponding to solutions of relativistic conformal fluid mechanics, were constructed 
in [12]. These solutions were obtained by starting from the general equilibrium configuration, 
promoting the parameters characterising it (temperature, relativistic velocities, etc.) to slowly 
varying functions of spacetime, and then solving the bulk field equations iteratively in a 
derivative expansion. A further non-relativistic limit leads to a correspondence between 
metrics of constant negative curvature and solutions of the incompressible non-relativistic 
Navier-Stokes equations of the underlying conformal fluid [13, 14]. 

In our construction the analogue of the thermal state is Rindler spacetime (see figure 1), 
and the first step in [2] was to obtain the general equilibrium configuration. Nearby solutions 
describing the approach to equilibrium were then obtained by promoting the parameters 
appearing in the equilibrium solution to slowing varying functions of spacetime. Finally, the 
bulk equations were solved iteratively by applying a non-relativistic hydrodynamic expansion. 
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which we wih call the e expansion. More precisely, the incompressible Navier-Stokes equations 
have a scaling symmetry (that, in particular, scales space and time non-relativistically) , and 
higher-derivative corrections to the Navier-Stokes equations are naturally organised according 
to their scaling. At third order in the e expansion, one finds a bulk metric corresponding to 
solutions of the incompressible non-relativistic Navier-Stokes equations [1], while at higher 
orders the bulk metric corresponds to solutions of the Navier-Stokes equations corrected by 
specific higher-derivative corrections [2]. 

In [2] it was observed that all information could be recovered from a relativistic dissipative 
fluid by taking a non-relativistic limit. As explained in [2], the e-expansion never decreases 
the number of derivatives, but it may increase it, so the complete answer up to a given order in 
e may be obtained starting from a relativistic dissipative stress tensor containing a sufficient 
number of dissipative terms but the converse is not in general true. Indeed, the relativistic 
expansion is considerably more compact: we found that almost the entire information up to 
order e' is encoded in just one first order and four second order coefficients, while two further 
second order transport coefficients were undetermined at this order in the e expansion. In 
fact, only two terms in the non-relativistic stress tensor up to order were not recovered, 
but both of these required starting from a relativistic stress tensor of third order in gradients. 

These results indicate that there is a manifestly relativistic construction, and one of the 
aims of this paper is to flesh this out: instead of solving the bulk equations iteratively in 
the non-relativistic hydrodynamic expansion, we will solve them in a relativistic derivative 
expansion. This new expansion is significantly more powerful, and allows a treatment of the 
entropy currents which would be hard without going to very high order in the previous non- 
relativistic expansion. We are also able to compute the previously undetermined coefficients 
in the second order expansion of the fluid stress tensor, as well as recovering our earlier results 
for the other coefficients. Thus the current treatment is the Ricci-flat analogue of the AdS 
treatment [12], while our previous construction is the analogue of [13, 14]. 

Let us now summarise the results for the stress tensor up to second order in gradients, 
which in the gauge where TabU^h'^ = takes the form 

Tab = PUaUb + phab + Ilfft, nff,u" = 0. (1.1) 

One of the main results is that the fluid dual to vacuum Einstein gravity has zero equilibrium 
energy density 

Peg = 0. (1.2) 

Moreover, the stress tensor, including dissipative terms, satisfies the quadratic constraint, 

dTabT"^ = T^. (1.3) 

This constraint determines p as a function of p and 11^^ and, as such, it may be considered 
as a generalised equation of state. When the relation is applied at equilibrium it leads to a 
quadratic equation with one of the two roots being (1.2). 
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The remaining freedom in defining the fluid variables is usually removed by redefining the 
energy density so that TahU°'vf' = p^q- This so-called Landau gauge cannot be reached here 
since the equilibrium energy density is zero. Instead, we take the isotropic gauge by imposing 
that 11;^^ does not contain terms proportional to hah- A general fluid in flat spacetime and 
at first order in gradients is determined by two first order coefficients: the shear viscosity 
and bulk viscosity. A short computation shows, however, that when the equilibrium energy 
density is zero the conservation of the stress tensor at leading order implies that the fiuid is 
incompressible (to this order). The bulk viscosity may then be replaced by another parameter 
C which measures variations of the energy density at first order in gradients, p^^^ = C,' D\xip. 
Up to second order in gradients, one needs 11 additional second order coefficients (in flat 
spacetime), namely 

p = C'D Inp + ^ (^dilCabJC"^ + d2^ab^"'^ + d:i{D \npf + d^DD Inp + d^iD^ Inp)^) , (1.4) 

nffe = -27^1Cab + ^ {ciKK^b + C2K.lS^\c\b) + C^^aVl^h + C^KhtdcOdiup 

+C5lCab D In p + ceD^ In pD^ In py (1.5) 

One of our main results in this paper is the computation of all the above coefficients for the 
fiuid dual to vacuum Einstein gravity, 

C = 0, di = -2, d2 = d3 = d4 = d5 = 0, 

r] = l, Ci = -2, C2 = C3 = C4 = C5 = -cq = -4 . (1.6) 

As noted earlier, for theories satisfying the constraint (1.3), the energy density (1.4) is not 
independent but rather follows from (1.3), so all in all the dual fiuid is determined by one 
first order (the shear viscosity) and six second order coefficients (the c coefficients). 

A general feature of systems away from equilibrium is that they possess an entropy cur- 
rent with non-negative divergence. At equilibrium this current should reduce to the conserved 
entropy current, = SeqU^, where Seq the entropy density at equilibrium. In the hydrody- 
namic regime, the entropy current may differ from this expression by terms of higher order 
in gradients. Here we classify the possible entropy currents with non-negative divergence for 
fiuids with vanishing equilibrium energy density up to second order in gradients (in fiat space- 
time) . It turns out that first order gradients are not allowed and that there is a five-parameter 
family of allowed entropy currents depending on second order gradients: 

J'^ = SeX(l + ^(ai/C,6/C"^ + a2nabn^' - ^(4a2 - 56i + 4^2 + 63) (^± Inp)')) 
+ ^ {hik'^^dblCl + b2DlD Inp + h^KlD^ Inp + (26i - 262 - hW^D^ Inp 

+(4o2 - 56i + 362 + b-i)Dl \npD Inp) . (1.7) 
A two-parameter subset of these entropy currents is in fact trivially conserved. 
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One may then ask what is the entropy current associated with the fluid dual to vac- 
uum Einstein gravity. To compute this we have adapted the construction of [30], defining 
the boundary entropy current by pulling back a suitable horizon quantity along certain null 
geodesies. This leads to an entropy current of the form (1.7), i.e., with non-negative diver- 
gence, with coefficients 

ai = 1, as = ^, h = -1, 62 = -2, 63 = 1- (1-8) 

Given the unconventional properties of the fluid dual to vacuum Einstein gravity, it is reas- 
suring that the entropy current indeed has non-negative divergence. 

It was observed in [1] that from the bulk perspective, the non-relativistic expansion could 
be expressed as a combination of a Weyl rescaling plus a particular near-horizon limit. A 
similar bulk interpretation also exists for the relativistic expansion we consider here. As one 
takes the near-horizon limit, however, there are different quantities that one keeps fixed in 
the two cases, so the two limits are distinct. Of course, as discussed earlier, one may always 
take a further non-relativistic limit to go from the relativistic to the non-relativistic case. 

This paper is organised as follows. In the next section, we present the relativistic con- 
struction of the near-equilibrium solutions. Then, in section 3, we present the solution to 
second order in gradients, and in section 4, we discuss the classification of entropy currents 
to second order in gradients and the holographic computation of the entropy current. In 
section 5, we discuss the near-horizon limits and we conclude in section 6. Finally, in the 
appendix we present a basis for scalars, vectors and tensors, up to the order required for the 
hydrodynamic analysis. 

Note added: During the completion of this paper we were informed about the forth- 
coming publication [31] which has significant overlap with the material presented here. The 
results of [31] are in agreement with those presented here. 

2 Relativistic construction of near-equilibrium solutions 
2.1 Seed metric 

We start with Minkowski spacetime in Rindler coordinates, 

ds^ = -rdr^ + 2dTdr + dxidx\ (2.1) 
The metric ^ab on the surface Sc defined by r = rc is 

-fabix'^dx^ = -rcdr^ + dxidx\ (2.2) 

where the coordinates x"" = (r, x*). To obtain the zeroth order seed metric we perform the 
following changes of coordinate: first, we send 

r ^ r + l/p^ - rc, r y/r^pr, (2.3) 
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taking the metric to 

ds^ = 2p{^cdT)dr -p\r- r,)(^dr)2 + -fabdx^dx^ (2.4) 
and second, we perform the boost 

y^r — )• — ltaX°, ^ — U^^/r^T + {1 + j)~^U^UjX-' (2-5) 

where ^ 

= f (-r-c, Vi), 7 = (1 - v^/rc)-'^^. (2.6) 



Vrc-v 

Since this boost preserves jabdx'^dx'', we arrive at the metric 

ds^ = -Ipuadx^dr + [jab - p\r - r^)uaUb]dx''dxK (2.7) 

A simple calculation reveals that, with this metric, the Brown- York stress tensor on Sc takes 
the form of a perfect fluid in equilibrium: 

Tab = Phab, Kb = lab + UaUb ■ (2.8) 

The position of the Rindler horizon is now ry^ = rc — l/p^, as may be seen by writing the 
metric in the form 

ds^ = -Ipuadx^'dr + [hab - p^{r - r^)nanb]dx''dx^ (2.9) 

In the remainder of this paper we will set Tc — ?• 1 (taking 7^6 — t- rjab)- This may be 
accomplished without loss of generality via the scaling 

(r, r, Xi,rc,p, Vi) (A^r, r, Xxi, X^rc, X~^p, Xvi). (2.10) 

Acting on the equilbrium metric this sends ds^ — t- A^ds^, after which the constant overall 
conformal factor may be dropped. One may similarly restore Tc at any point by reversing 
this procedure. 

2.2 Integration scheme 

We start with the zeroth order seed metric (2.7) in the form 

ds^ = -2puadx''dr + [r]ab + (1 - e)uaUb]dx''dx\ 6* = 1 + ^^(r - 1) (2.11) 

where we have scaled rc to unity and introduced the quantity 6 for our later convenience. The 
position of the horizon is now r-^ = 1 — 1/p^. Note also that the relativistic fluid velocity is 
normalised such that rj'^^UaUi, = — 1. In this metric, the velocity and pressure Ua and p should 
now be regarded as functions of x'^ = (r, x*). The inverse of this metric is 

5- = r-rH, 5™ = -n^ ff'^" = (2.12) 

P 
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where we define 

f^ab _ ^ab ^ ^a^b^ _ ^ab^^_ ^2.13) 

Weighting derivatives such that dr ~ 1 and da ~ e, adding a piece q^m) to the metric at 
order engenders a change in the Ricci tensor 

1 

—I 
2 



2p 



1 



2p 
1 

-^UaUb 



p'^ ^ ' p 



(2.14) 



A convenient gauge choice is g)-^ = for n > 1. This choice eliminates metric components 
for which we do not have natural boundary conditions, and moreover ensures that worldlines 
of constant are bulk null geodesic to all order (see section 4.2). The linearised Ricci tensor 
tensor is then 



2''' "riJab ' 



- J_ 1,6(52 Jn) P jhcg in) 
"^ra ~ 2p ^^"-^ ^'J'a'l- 'Jrghr. i 



(2.15) 



(n) ~ (n) 

Setting = 6Rliu + RliJ , we obtain the integrability conditions 

= dr l^h'^'R^^' - ^ Rl^^) - R\^^ > = OR^-J + pv'R^:^ . (2.16) 

The first of these is the r-component of the Bianchi identity at order e^] while for the second, 
the a-component of the Bianchi identity enforces 



Evaluating the Gauss-Codazzi identity on Sc at order e*^, we find 



(2.17) 



V'T.^") = [2V\K^ab - Kab)t'^ = [-2i?,^iV^](") = —fi^\x). (2.18) 

Thus, conservation of the Brown- York stress tensor at each order ensures that one may 
integrate the bulk equations. 

The radial Hamiltonian constraint evaluated on Sc is given by 



- KabK'^' = 0, 



(2.19) 
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which upon using the definition of the Brown- York stress tensor, Tab = %^{^lab — Kah), 
becomes the quadratic constraint given earher in (1.3). As was discussed in [2] and reviewed 
in the introduction, this constraint plays the role of a generalised equation of state: given p 
and lij;^^ it determines p. 

Note that the radial Hamiltonian constraint in AdS plays exactly the same role, i.e., it 
determines the equation of state. Indeed, in the presence of a cosmological constant the radial 
Hamiltonian constraint (in Fefferman-Graham gauge) becomes 

K"^ - KabK"^ = d{d - 1), (2.20) 

where have set the AdS radius to unity and assumed a flat boundary metric. Expanding 
about conformal infinity, the extrinsic curvature is given by [32] 

Kab = riab + i^(d)a6 + • • • , (2-21) 

where the subscript indicates dilatation weight and the dots represent higher order terms 
that do not contribute when we evaluate the constraint at conformal infinity. Inserting the 
holographic stress tensor [32] 

Tab = 2{K(^d)riab - K(^d)ab), (2.22) 

into the AdS Hamiltonian constraint (2.20) then yields 

T^ = 0, (2.23) 

which implies the equation of state of a conformal fluid. ^ 

Returning to the case of the Rindler fluid, a particular integral of (2.15) is 

~g^^ = a^-)uaub + 2/3f:)^,) + 7!?, u^~P^^) = 0, ^^^7!? = 0, (2.24) 



where 



aW = ci")(x) + (1 - r)c^^\x) + 2p^ C dr' C Ar"{h^''R^^^ - ^R^^^ 



2 



= ctJix) + (1 - r)c4a(x) + 2p [ dr' [ dr"/.^i?(;!\ 

Jr Jr' 

7i? = 4:l(x) + (rr) In^ - 2/ dr'i f' dr"/^>fi?(?, (2.25) 
where = g^^'^^'" R^^J . Note that to satisfy the rr equation, we must have that 

2^W = h^'d^,^^^ = -P^h^'c^:J, + - ^ £ dr'h^'R^:;^. (2.26) 

^ If we instead consider a general boundary metric then the r.h.s. of (2.23) contains the holographic Weyl 
anomaly. 
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Using the first integrability condition in the form 



(2.27) 



one can show that 



eh'^'R^:^ - ^ + \h^'R^:i - ^ d//.^^i?i?, (2.28) 



where [...]* indicates evaluating at r = r*. One then obtains an equation for r*, namely 



abin) 



6ab p2 



(2.29) 



Later, we will choose to set the coefficients of all logarithmic terms to zero. For the trace 
component above, this entails setting the lower limit of integration such that = r-^ = 
1 — whereupon 6^, = 0. 

Allowing now for gauge transformations ^("^^ at order e", as well as re-definitions (5u^"^"(x) 
and {x) of the fluid velocity and pressure, the solution above generalises to 



gil^ = -naipdr^^"^ - (1 - 9)ntdr&^' + ,5p(")] + r^abdr&^' - p5u^:'\ 



in) 

b) ■ 



9ab = ~€b - uau,[p\^^> + 2p5p^^\r - 1)] + 2(1 - e)u(M, 
To impose the gauge choice gi'j} = 0, we must then set 

p 

where Ua^jf^"" = and 

^(n) ^ |(n)(^)^ ^(n)a = _(l _ r)p5u^^^'' + ^^"(x). 

The remaining metric components then take the form: 

aW = cj") ^ ^1 _ ^^(^C«) + p^p{n)^ + 2/ d/ dr"(/i^'^i2^;j) - -^(")) 

P'-a^ = 4? + (1 - 0(c4a +p'<5ni")) + 2p dr' dr"h\R^ 

J r J r' 



(2.30) 
(2.31) 
(2.32) 



(2.33) 



(2.34) 



„{") _ J") I J^l 



€Jb + cZl^0-2l dr' 



r - r-H 

(n) 



dr"KhtR^-\ 



(2.35) 



Imposing the boundary condition g^^^^' = for n > 1 on X]^ then fixes 



C3a' = 0, = 0. 



(2.36) 
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Moreover, for regularity on the future horizon r = r-^, we must set 

4al = 0, n = rn. (2.37) 

In summary then, our integration scheme is the following. The rr and ra metric compo- 
nents are given to all orders by 

9rr = 0, g-ra = -pUa, (2-38) 

while the ab metric components may be decomposed in fluid variables as 

,i^)=aW^,n, + 2/3f:)n,)+7l?. (2.39) 

Beginning with the seed metric g^jj = —6uaUh + hab, the metric at all subsequent orders is 
then given by 

= (1 - r)F(")(x) + 2p2 [\r' C dr" {h^'^ R^^^ - ^i?^")), 

Jr Jr' 2 

= (1 - r)F!r\x) + 2p f dr' dr'n'^R^^ , 

Jr Jr' 

7i? = -2 dr'-l- r' Ar"KhiE!^I^l (2.40) 



r' - rn 



Here, the arbitrary functions 

F(")(x) = 4")(2;) + p8p^''\x), Fi")(x) = C4a(x) + p^8u^^\x), (2.41) 
encode the choice of gauge for the dual fluid, and will be fixed as we discuss in the following 

(n) (n) 

section. Note also that Fa is transverse: u°'Fa = 0. 
2.3 The Brovi^n-York stress tensor 

(n) 

The variation in the extrinsic curvature of Sc at order e"" due to g^^ is 

= \£Ng':^ = \N^dj:^ = (2.42) 

(Note here that the normal A^'^Iec = P~^8r + 5a it" to all orders, since the bulk metric at Xl^ 
is effectively fixed.) Evaluating this explicitly, we find 

^""^^k = -^^^'^^^'^-^ - ^C^^^ dr'/.>^i?(?. (2.43) 

(n) 

The variation in the Brown- York stress tensor due to is thus 

<^rl"\^ = 2(.?„,5i^(")-<^K(^)) 

= -F^^^hab + -HaF^:^ + 2p [ dr'i^Mh''' - Kht)Rt^. (2.44) 
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The complete Br own- York stress tensor on Sc at order e" is then 

T:;^^k = ^T^^+T^\ (2-45) 
where T^^"* represents the contribution at order e" due to the metric up to order e""^. 
2.4 Gauge choices for fluid 

We will define the relativistic fluid velocity u"' such that 

= (2.46) 

(n) 

which then uniquely fixes Fa ■ 

Fjr^=phifl:\^. (2.47) 

To fix F^^\ we impose that there are no corrections to the pressure, i.e., that the coeffi- 
cient of hab in Tab is fixed to be exactly p. From (2.44), we see that F^") is then determined 
uniquely. 

3 Solution 

In the previous section, we saw how to systematically construct the near-equilibrium solution 
in terms of a relativistic gradient expansion starting from the seed solution. We saw moreover 
that the solution is unique once we impose the bulk gauge conditions, the gauge conditions on 
the fluid stress tensor, and regularity at each order in the expansion. In the present section, 
we will now explicitly compute this solution to second order. 

3.1 First order 

Computing the Ricci curvature of the seed metric, we obtain 

= 0, i?W = 0, b!-^^ = {Dp + ^d,u')uaUb + pu^aab) + U^adb)P, (3.1) 

where D = u^da, = h\db and the acceleration Oc = Duc- Since h'^^R^b vanishes, the 
integration step is trivial and we have 

gl^,^ = {l-r)[F('^UaUb + 2Fiyb)]. (3.2) 
Evaluating the Brown- York stress tensor on Sc, we flnd 

Tab\^^ = {p + 2D Inp + ^f(i)) hab + 2u^a {^ab) + ^Flf") - 2ICab, (3.3) 
where we write the fluid shear and vorticity 

K^ab = Kh'^d(^^Ud), Qab = Khfd^^Ud], (3.4) 
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so that 

daUb = ICab + ^ab " UaO-b- (3.5) 

Note here that conservation of the stress tensor at zeroth order yields the conditions 

dcu" = 0{d'^), ac + D^\np = 0{d^), (3.6) 

and that we used both of these to simphfy the form of the stress tensor at first order. 
Our gauge condition for the pressure immediately sets 

F(^) = -2Dp, (3.7) 
while that for the fluid velocity at first order reads 

= /i^TbX = -^fP - 2aa F^i) = -2paa. (3.8) 

In conclusion then, 

9ab = "^(^ ~ ^){uaUbDp + 2pa(^aUb)) (3.9) 

while 

Tab\j^^ = phab - 2ICab + 0(3^), (3.10) 

from which we may read off the viscosity rj = 1. We also deduce that the first order coefficient 
= 0, where C' is defined in (1.1)-(1.4). 

3.2 Second order 

The fluid equations of motion at second order may be directly obtained from conservation of 
the first order stress tensor (1.1) with ^' = 0, ?7 = 1 as 

dbu^ --KbcKf"" = Oid^) , aa + Djr\iip--hldbK}' = Oid^). (3.11) 
P P 

Next, from the first order metric (3.9), one may obtain the second order piece of the Ricci 
tensor R]iy . The computation is straightforward but laborious. Using (A. 3) - (A. 10) and 
the fiuid equations (3.11), the result may be expressed in the basis of tensors given in the 
appendix. We find 

R?} = PUa ({D^ Inpf + ICbcIC'' +p\r- mbc^'') + pKdblC\ + p{ICab + J^afc)!?! Inp , 



R^ab = UaUbO{r) [{D^ \npf + tCcdlC"" + p\r - + 2u^ J{r) ((/C^), + Ofe),)!?! Inp 

+ hl^ddJC"^))^ + 2JCabD Inp + 2hlhidcdd Inp - 2D^ lupD^ Inp + JC^IC^b + K'n,b 
+ 2/C,(,f]^,) + 2p\r - m^ia^^) ■ (3.12) 
(Note the equations N^'R^a = 0(9^) are satisfied, as required by (2.18).) 
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(2) 

As noted earlier, the stress tensor at second order T^^ evaluated on Sc is the sum of two 

(2) (2) 

pieces: the piece ST^^^ resulting from the second order contribution to the metric g^^^ , and 

" (2) 

the piece T^^ resulting from evaluating the stress tensor to second order using the metric up 
to first order. Computing the stress tensor at first order without using the first order fiuid 
equations of motion, we obtain 

Tab = Phab - + 2)C{hab " UaUb) - 2n(„ (ofc) + D^^ lup) + 0(8^) . (3.13) 

Using the equations of motion at second order (3.11), we deduce that 

t!^^ = ^'^cdlC'^'iKb - UaU,) - ^u^, {hlddlC'')^ . (3.14) 

Next, inserting the second order Ricci tensor (3.12) into (2.44), and making use of the sub- 
stitutions listed in the appendix, we find 



^T^b = Kb {^-F^^^ - ^(^ab/C-')) + ^n,nb(/C,6/C" 

■ (hf^Di^Di \np - \npDi \np - KhtDl^a + 2/C,(,J]'=,) - 1^,(^0^,)) . (3.15) 



2 

P 

Combining these two contributions to the second order stress tensor, our fiuid gauge conditions 
in section 2.4 are met when 

= -2/Ca6/C"^ = 2hlddlCi . (3.16) 

The fluid stress tensor then takes the expected form (1.1), with p and 11^^ as given in (1.4)- 
(1.5) and with the advertised second order coefficients (1.6). 
As a consistency check, we note that our expression for p^^^ , 

p(^) = --JCabJC""" , (3.17) 
P 

as well as the coefficients ci, C2, C3 and C4 listed in (1.6) coincide with our earlier results 
obtained in [2] using the non-relativistic e-expansion to order O(e^). The two remaining 
coefficients, C5 and cg, that were not determined in the analysis of [2] have moreover now 
been computed. 

Recalling that the second order metric takes the form 

g^^ = a^'^^UaUb + 2f3[y,^ + 7^,^ , (3.18) 

using the solution algorithm (2.40) with our result (3.12) as input, we obtain 

= 2(r - l)/C„,/C'^^ +^2(1/^,/^^^ + {D^ lnpf){r -lf + ^{r- ifQab^''' , 
/3(2) = -2(r - l)hlddlC'', + p\r - l)\hld,IC% + (ICab + f^ab)^l Inp) , 
7^,) = 2(r - 1) (2hlhtdcdd Inp + 2JCabD Inp - 2D^ InpD^ Inp + JC^IC^b 

+ 2n,'na + 2/Ce(af^"b)) + p\r - ifn.an^h . (3.19) 
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As a check of this result, one may proceed to expand the metric we have just calculated to 
O(e^) in the non-relativistic e-expansion. One should then recover all terms of up to second 
order in gradients in the metric of [2], modulo a small complication which is that the metric 
in [2] was computed in a slightly different radial gauge. (The choice of fluid gauge in [2] is 
also slightly different.) To get round this, we simply repeated the analysis of [2] with our 

(2) 

current gauge choices. Comparing with the e-expansion of , we then found that the two 
results are indeed consistent. 

In summary, the complete metric up to second order is given by 



9r 



where 



0, 



a 

Pa 
lab 



9r 



-pUa 



gab = aUaUb + 2/3(aUf,) + jab , 



(3.20) 



— 1 — p^(r — 1) + 2pD lnp(r — 1) + a 
+ 2pa,(r-l) + /3(2), 



(2) 



(3.21) 



The inverse metric up to second order is 



1 



5 = 



a + h'^'Mb) 



(3.22) 



3.3 Fluid divergence at third order 

Given the stress tensor at second order, it is straightforward to obtain the fluid equations of 
motion at third order. We will need to make use of the third order fluid continuity equation 
u°'d^Tab = in our forthcoming discussion of the entropy current. Let us then derive this 
equation for a general fluid stress tensor of the form (1.1) in combination with (1.4)-(1.5). 
Using the basis of tensors given in the appendix, we obtain 



daU'' = - 
P 



1 r 
+ — 



C4 + 2di)IC''''dadblnp+{-ci + 2di)/C"''/C;j/C 



abt 



be 



+(-C3 + 2di - 4:d2)IC'"'na''nbc + 3d2nab^l'"'D\np + (-C6 - 2di)/C"''L>^ InpD^ Inp 



+(C4 - C5 + di - 2r]C)D InplCablC'' - d^D^ \np + {-2d^ + di + (^)D^ InpD Inp 



+d3{D\np)'^ + d^iDjL lnp)'^Dlnp - 2d^DD\^ Inpi^f Inp + 0((9^) . 
Interestingly, upon inserting the coefficients (1.6) this equation simplifies to 



(3.23) 



daU'' = -JCabJC"' - ^{IC^'lCalCbc + I? Inp/Cfe/C"") + 0(8*) 

p 

2 



P 



"ab^ 



1 

-I 

p 



Dlnp . 

p 



+ o(a^ 



(3.24) 
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In particular, the right-hand side is non-negative since the third order terms are small cor- 
rections to the non-negative second order term. This is a special property of the fluid dual 
to vacuum gravity: in the general case, the divergence might be negative for flows where ICab 
is small and other terms in the expression (3.23) dominate. 

4 Entropy current 

4.1 General fluid entropy current 

A general entropy current is by deflnition a current which has non-negative divergence given 
the fluid equations of motion. We will derive in this section the constraints on the form of 
the entropy current at first and second order in the derivative expansion using solely the fluid 
equations of motion (3.11). In particular, our considerations do not depend on the second 
order part of the bulk metric. 

Using dimensional analysis and the classification of first order scalars and vectors given 
in the appendix, the general entropy current at first order has the form 

J'' = Seqiu'' + -u''D\np+^Dl\ni^ +0(9^). (4.1) 
\ P P J 

Here, we normalised the current using the equilibrium entropy density s^q = 1/(4G). Using 
the relationships in the appendix, the divergence of this entropy current is 

—daJ" = '^ICablC'^' + ^natn""' + -{D^ Inp - [Dlupf) - ^{D^ \npf + 0{d^). (4.2) 

Seq P P P P 

Since the two terms Q.ah^"'^ > and — (D^lnp)^ < of opposite sign might dominate the 
divergence when the shear tensor is small, we need /3 = 0. Also, since the term D'^lnp — 
(Dlnp)^ of indefinite sign might dominate the divergence when the shear tensor is small, we 
need a = 0. The entropy current J" = Sequ"" is therefore the only possible expression at first 
order. 

From the classification of second order scalars and vectors in the appendix, the general 
entropy current at second order takes the form 

J'' = SeX(l + ^{ailCablC^ + asl^ab^^"^ + as{Dlnp)^ + a^oHiip + a^iD^ Inp)^)) 

+ ^ [bih^^dblCl + b2DlD Inp + bslC^D'i Inp + hiVL^^Di Inp + h^Dl InpD Inp) , (4.3) 

where the coefficients a,, 6j are constrained by the fact that the divergence of the current is 
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non-negative. This divergence is given by 

daT = —JCatJC"^ + ^ (iC^D^Dt lnp(-2ai - h + 462 + 63) + K,tK!'JCl{-2 - 2ai + 262) 
p V 

+ ICti^\i}'^i-2ai + 4a2 - ih + 662) + I? Inpfi„fe0"''(-4a2 + 36i - 62 + 64 + 65) 
+ dalC^Dl lnp{-2bi + 262 + 63 + &4) + (I? Inp)3(-2a3) 
+ D In plCablC'' {-2 - 2ai - 62 - ^5) + D'' lnp(a4) + D'^ InpD lnp(2a3 - 204) 
+ IC'^D^ \npDi lnp(2ai + 61 - 862 - 263 + &4 + ^5) 

+ {D^ InpfD lnp(-2a5 - 62 - 65) + ^^^1 lnpZ?f lnp(2a5 + 62 + 65)) + 0(5^)- (4.4) 

The first term is leading in gradients and therefore the current has non-negative divergence 
in most cases. It might turn out, however, that the shear tensor is small and then third order 
gradients might be leading. Of these third order terms, some might form a perfect square 
with the term ^/Cafe^"^ when higher order gradients are taken into account, leading to a 
non-negative contribution to the divergence. There are other terms, however, which clearly 
cannot form a perfect square and which may be large even though the shear is small. These 
latter terms must therefore vanish giving rise to the necessary positivity conditions 



= -4a2 + 36i -62 + 64 + 65, (4.5) 

= -26i + 262 + 63 + 64, (4.6) 

= a3, (4.7) 

= a4, (4.8) 

= 205 + 62 + 65. (4.9) 



Further conditions might be found by studying the constraints up to third order or by deriving 
the constraints associated with putting the fluid in a curved spacetime. We will not perform 
such an analysis here. The constraints (4.5)-(4.9) should be obeyed by any physical entropy 
current. Using these constraints to eliminate 03, 04, 05, 64 and 65 in (4.3) leads immediately 
to the five-parameter family of entropy currents with non-negative divergence (1.7) given in 
the introduction. 

As a final remark, we note that certain entropy currents may be written as the divergence 
of an anti-symmetric potential, i.e., 

J« = dbX^"''^ . (4.10) 

Entropy currents of this form are trivially conserved and describe lower-dimensional con- 
servation laws at the boundary of the fluid domain. Since we consider only an infinitely 
extended domain, we will ignore such boundary terms. The trivial entropy currents may then 
be straightforwardly classified as 

J- = db (ti'-^n'^' + ts^nl'^D? Inp) . (4.11) 
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Making use of the relations in the appendix, one finds that the trivial entropy currents take 
the general form (4.3) with 

oi = —t2, 02 = —h + 12, as = 04 = 0, 05 = -t2, 
bi = -ti, b2 = -t2, b3 = -ti + 2t2, bi = -h, 65 = 3*2, (4.12) 

which obviously satisfy the constraints (4.5)-(4.9). Removing these trivial entropy currents, 
one can reduce the five-parameter family of entropy currents (1.7) to a three-parameter family 
of non-trivial entropy currents. 

4.2 Defining the holographic entropy current 

Our bulk metric takes the form 

ds^ = -2p{x)ua{x)dx''dr + gab{r, x)d3;"dx^ (4.13) 

implying wordlines of constant are null geodesies to all orders. These null geodesies define 
a natural map between points on the horizon and on the boundary. In the spirit of [30], we 
will obtain a boundary entropy current by pulling back a suitable horizon quantity along these 
null geodesies. In this section we present our prescription in detail and derive a formula for the 
entropy current which extends that of [33, 34] to allow the use of a non-affinely parametrised 
horizon generator. 

We begin by foliating the bulk spacetime with a family of null hypersurfaces defined by 
S{r,x) = constant, with the horizon corresponding to the surface S{r-}{{x),x) = 0. In the 
equilibrium case, the horizon is defined to be the Rindler horizon of the bulk solution. In the 
near-equilibrium case, we will assume that the fluid is returned to equilibrium in the limit 
of late times through the action of dissipative forces. The horizon of the near-equilibrium 
solution is then defined as the unique null hypersurface which asymptotes to the Rindler 
horizon of the late-time equilibrium solution. Expanding the location of the horizon r-^(x) in 
fluid gradients, the zeroth order term must therefore match the location of the Rindler horizon 
of the equilibrium solution, giving ry_{x) = 1 — 1/p^ -|- 0{d). The higher order corrections 
may then be obtained by requiring the horizon to be null; we will return to solve for these in 
the next subsection. 

On the horizon then, we have 

^=^a\n=[^rSda'^ + daS]^. (4.14) 

An afhnely parametrised normal vector to our family of null hypersurfaces is £^ = d^S, since 

^"^u^^. = \d^{^') = 0. (4.15) 
The vanishing of i"^ = i'^df^S everywhere also implies 



r = -{drSrH^'daS. 



(4.16) 



Denoting the expansion evaluated on the horizon by ^-^ = (V^/'^)-^, we have 

= daiiV^nn] - V^n^ndr{idrSr'daS)\n, (4.17) 

where to obtain the second hne we used the chain rule in the form da[{\/—gi"')n] = 9a{\^ — g(^"')\n + 

driV^nndarH- 

Close to the horizon, we may Taylor expand S{r, x) as 



5(r, x) = {r- rnix))Si{x) + -(r - rn{x)fS2{x) + 0(r - rnf, 



(4.18) 



where the absence of a zeroth order term is required by 5(r-^,x) = 0. A simple calculation 
shows that Si is nonzero for the Rindler horizon of the equilibrium solution; since corrections 
at higher order in gradients cannot cancel this leading term, Si is everywhere nonzero in the 
near-equilibrium case as well. In terms of this Taylor expansion, we obtain the exact relations 



where 



dr{{drS)-^daS)\n = da hi^i, = Siig''' - g^^Obrnjn 



e = 9^''d,{T-rH{x)) 



(4.19) 



(4.20) 



and we have chosen l':^ to be future-directed ensuring that Si > 0. (Note is future-directed 
for the equilibrium solution, and hence for the near-equilibrium solution also.) Since £^ = 0, 
we must have = 0. Solving this latter condition in the hydrodynamic gradient expansion 
provides us with the location of the horizon r-^(x), to which we will return shortly. 
Combining (4.17) and (4.19), we have 



V^n^n = daiV^n^^) " V^n^udalnSi = Sida 

and hence, introducing the boundary metric g^, we may write 

V^n^H 1 



Si 



where the entropy current 



-5e 



5a(^/^s^")=vPJ^ 



ja 



1 



~9h 



Prom the Raychaudhuri equation, 

1 



4G7V ^/^s 



d 



n - crabO'"'^\n < 0, 



(4.21) 



(4.22) 



(4.23) 



(4.24) 
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where the dot denotes the derivative with respect to the affine parameter along the horizon 
and aab is the shear of the geodesic congruence, while the vorticity vanishes since is 
hypersurface orthogonal. Since the fluid returns to equilibrium in the limit of late times, and 
0-^ = for the equilibrium solution, it therefore follows that in the near-equilibrium case 
^"H ^ at all times. (Note this conclusion relies on the fact that Oy, is the expansion with 
respect to the affine generator iij^. If instead one tried to use the expansion defined with 
respect to the non-affine generator then the Raychaudhuri equation would acquire extra 
terms of indefinite sign, invalidating the argument.) 

Given then that 0-^ is non-negative and Si > 0, from (4.22) the divergence of the entropy 
current must also be non-negative: 

y(^S)^a > Q_ ^4 25) 

Examining (4.23), we note that while the expansion O-y is necessarily that of the affinely 
parametrised generator if^, the current ^7*^ may nevertheless be expressed in terms of the 
non-affinely parametrised generator S^y^. 

Finally, let us discuss briefly two potential sources of ambiguity in the definition of the 
holographic entropy current. Firstly, pulling back to the boundary along a different set of bulk 
null geodesies will lead to a different boundary entropy current. Such ambiguities have been 
discussed in [30] and correspond boundary to boundary diffeomorphisms. A second potential 
source of ambiguity, mentioned in [30] and discussed in [33, 34], concerns the choice of bulk 
horizon (such as apparent horizon, etc.). We leave further investigation of these interesting 
issues to future work, and in the following, we focus exclusively on the entropy current defined 
in (4.23). 

4.3 Location of the horizon 

To evaluate the entropy current according to our formula (4.23), the first step is to compute 
r-^(x), the location of the horizon. As explained above, this follows from solving the null 
condition = 0, which from (4.20) reads 

= /^(r«) - 2g''^{rn)darn + 9''\rn)darnd,rn. (4.26) 
The solution takes the form of a gradient expansion 

r^(x)=r5J)+r«+rg) + ... (4.27) 
where = 1 — 1 /p'^ is the equilibrium position of the horizon obtained by solving 

5(o)(^l?^)=P^' + ^S^-l = 0, (4.28) 

and rj^^ contains terms of n-th order in gradients. At each order n in gradients, the equation 
reduces to a linear problem due to the fact that the only term involving is g'^^^ {vy^^ ) = . 
At first order, (4.26) reads 

= r« - ^a(^)(r2^) - -u^dj^^ , (4.29) 
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with solution 



n-^D\np. (4.30) 



At second order, we obtain the equation 



- ^Dl InpD^rfl^ + h'^'' davf^ d^^rf^^ , (4.31) 



with solution 



rJJ) = 1 (^ADD Inp- 8{D Inpf - ^JCabJC'^' - ^l^a6^^"' + {D^ Inp)' j . (4.32) 

Having obtained r-}{{x), we may now evaluate the non-affine horizon generator ac- 
cording to (4.20), making use of the inverse metric (3.22) evaluated on the horion. We find 

= ^ + ^(2DllnpDlnp - 2DlDlnp - h'^'^dcJCl + {IC\ + n\)D''^lnp^ + Oid"^). (4.33) 

We will not need S^y^ in what follows. The piece of ^-^ normal to the fluid velocity is simply 
the equilibrium term, with no corrections at either first or second order. The remaining piece 
of S^y^ tangent to the fluid velocity has only second order corrections. 

4.4 Evaluating the holographic entropy current 

In our case the boundary metric is simply Minkowski and l/4GAr = 47r, so the entropy current 
(4.23) reduces to 

J- = 4vr^^e^. (4.34) 

With ^y given in (4.33) above, it remains only to evaluate the determinant factor ^J—g-y- 
As an initial step, we first evaluate the determinant of the seed metric 

dsjo) = -2puadx''<lr + QabAx^dx^, (4.35) 

where, in the above and in the following, for clarity we will temporarily write 

Qab = gfb = "Hab -p^{r - l)UaUb . (4.36) 

Since 



—pu^ \ I 1 —pu^ \ I — p^u^g -"^u 
-pu g \ g M -pg^^u 1 



we find 

-ab 



(4.37) 

det 5(0) = -p'uaUbg"" det g, (4.38) 
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where the inverse metric 

(1 +p^(r - 1)) 

and the determinant det ^ = —(1 +p^{r — 1)) may be evaluated by the usual formula. The 
seed metric therefore has determinant 

det5(o) = V- (4.40) 

The determinant of the full metric up to second order in gradients may now be obtained 
perturbatively. Writing 

9iiu = g{o)^,u + 9{i)f,u + 9(2)iiu + 0{d^) (4.41) 
and expanding the formula det g = exp(tr In g), we find 

detff = det 5(0) (^l + ^(^^(i)/.!/ + 9^o)9i2)^lu - ^9'^o)9ii)up9lo)9{i)af, + ^i9^o)9ii)f,uf^ + 0(<9^). 

(4.42) 

Since in addition, 

9{i)rtJ. = 9{2)rtJ. = 0, 5(0) = h"-'', 7(l)a6 = Khh9{l)cd = 0, (4.43) 
making use of (4.40), we find simply 

det5 = V(l + h'^\2)ab) + 0{d^). (4.44) 
Evaluating this formula on the horizon r-^ = 1 — + 0{d), we obtain 

V^-H =P+ -iCabic"' + ^n^bn'"'' + o{d^), (4.45) 

p Zp 



and thus the entropy current 

J- = attu'^ (i + l/Cfce/c^^ + ^^fibcn'^ 

\ p'^ Ip^ 

47r 



+ ^ [2DI InpDlnp- 2DlD\np - h^^'dcIC^ + m + n"'^)D'i Inpj + 0{d^). (4.46) 
The entropy current takes the general form (4.3) with coefficients 

ai = 1, ^2 = ^, 03 = 04 = 05 = 0, -61 = 63 = 64 = 1, -62 = 65 = 2. (4.47) 



From (4.4), we obtain 

L 

P p^\ 



+ SIC^^'D^ \npD^ Inp - ADlnplCablC'') +0{d^). (4.48) 



All five entropy conditions (4.5)-(4.9) are obeyed, confirming that the divergence of the en- 
tropy current is non-negative as expected. 
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5 Near-horizon limits 



It was observed in [1] that the non-relativistic hydrodynamic expansion can be expressed as a 
near-horizon hmit when combined with a specific Weyl rescaling. In this section we show that 
the relativistic expansion can also be expressed as an alternative near-horizon limit combined 
with a Weyl rescaling. 

We noticed in section 2.1 that the equilibrium solution admits the scaling transformation 
(2.10) that generates a global Weyl rescaling of the metric. This Weyl rescaling still exists for 
the complete metric with higher-derivative corrections that we found in section 3 (with the 
factors of restored). Indeed, the metric has coordinates r,T,Xi and parameters rc, p{t,x) 
and Vi{T,x), where the relativistic velocity is decomposed as in (2.6). Equivalently, one 
can express the solution using the position of the horizon r-^(r, x) instead of the pressure 
after inverting the relation (4.27). The scaling 

{r,T,Xi,rc,rn,Vi) (A^ r, r, A Xj, rc, r-^, A Wj) (5.1) 

is equivalent to a Weyl rescaling 

ds^^A^ds^ (5.2) 

of the full near-equilibrium metric. Since we are interested in Ricci-flat metrics this constant 
overall factor may be dropped. 

We now want to consider the near-horizon limit Tc — t- r-^ — t- while preserving r-^ < rc 
and Tc — > (recall that ^/r^ also plays the role of the speed of light). Thus, in general 
we must scale the parameters 

{rc,rn,Vi)^ {Xrc,~X''rn,~X^v^) (5.3) 

such that a > 1 and b > 1/2. According to this point of view, any suitable choice of a and b 
defines an acceptable near- horizon limit. 

Recall that the non-relativistic e-expansion - defined as the homogenous scaling trans- 
formation of the incompressible Navier-Stokes equations - is given by 

{r,T,Xi,rc,rn,Vi) ^ (^r, ^,^,rc,e^rn,evi^ (5.4) 

Recall also that the incompressible Navier-Stokes equations are an attractor under the e 
scaling in the sense that when e — t- all higher order corrections become small. 
Combining this expansion with a Weyl rescaling (5.1), setting A = e, leads to 

{r,T,Xi,rc,r-u,Vi) ^ (^e^r, ^, x^, e^r^, e^r-^, e^Ui^ , (5.5) 

which defines a near-horizon limit (5.3) with A = and a = 2, 6 = 1, i.e., we consider the 
limit 

Tc — )■ 0, = fixed, — = fixed. (5.6) 
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Let us now consider the relativistic limit 

{r,T,Xi,rc,rn,Vi) (r, j, ^,rc,ry,,v^ , (5.7) 

i.e., the pressure p and u'^ are zeroth order quantities, only derivatives carry weight ~ e. 
The ideal relativistic fluid equations (3.6) are an attractor under the relativistic scaling in 
the sense that when e — )• 0, all higher order corrections become small. 

Combining this expansion with a Weyl rescaling (5.1) with A = e leads to the near-horizon 
limit 

{r,T,Xi,rc,rn,Vi) (^e'^r,j,Xi,frc,e'^rn,m^ , (5.8) 
which is (5.3) with A = and a = 1 and b = 1/2, i.e., we consider the limit 

2 

Tc — )• 0, — = fixed, — = fixed. (5.9) 

In particular, this means that we keep fixed relativistic velocities as yjrl plays the role of the 
speed of light. Notice that under (5.8) the normalised horizon generator C, = , ^ ^ {dr + Vidi) 
(where C^lsc — see [2]), is invariant and the temperature and pressure satisfy a simple 
scaling law T — )■ T/e and p ^ p/e, while the transformations under the non-relativistic scaling 
(5.5) are more complicated. 

6 Conclusions 

In this paper we presented a construction of a (d + 2)-dimensional Ricci-flat metric corre- 
sponding to a (d + l)-dimensional relativistic fluid with specific transport coefficients. In a 
specific non-relativistic limit we recover the results discussed in our previous work [2]. We 
have further obtained a holographic entropy current with a non-negative divergence, in ac- 
cordance with the second law of thermodynamics. We also showed how to reinterpret the 
relativistic hydrodynamic expansion as certain near- horizon limit. 

There are many interesting directions that one may wish to pursue further. Some of the 
numerous questions that were raised in [2] have now been addressed, both in the literature 
discussed in the introduction, and in the present work. Nonetheless, many interesting ques- 
tions remain. Perhaps most far-reaching of these are the questions concerning holography. 
How concrete can we make this holographic duality? Can we move away from the hydro- 
dynamic regime? Can we set up holography for general spacetimes by using the discussion 
here as a local holographic reconstruction of small neighbourhoods which should then be 
patched together to obtain a global description? Answering any of these questions would be 
a significant step towards formulating a general theory of holography. 
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A Basis of hydrodynamic scalars and vectors 

At equilibrium, the long wavelength dynamics of relativistic quantum field theories in flat 
spacetime at vanishing charge density can usually be described by a set of fundamental fluid 
variables consisting of the energy density p, pressure p and the fluid vector u"", constrained 
by an equation of state p = p{p). For many applications it is usually more convenient and 
natural to trade the pressure for the temperature T and treat (T, u"") as the fundamental 
variables. At non-zero energy density either choice of fundamental variables is equally valid 
but the equation of state of the equilibrium fluid dual to vacuum Einstein gravity is p = 0. 
Conservation of the stress tensor Tab = phah + 0{d) then leads to the incompressible ideal 
fluid equations 

dau"" = ^ + 0{d'^), a" = -D^\np + 0{d'^). (A.l) 

The fluid is incompressible at flrst order in gradients precisely because the equilibrium energy 
density is zero. As a consequence, the bases of hydrodynamic scalars and vectors tradition- 
ally used to describe the fluid dynamics at higher orders in gradients, see e.g. [35], are not 
applicable to this special case. In this appendix we construct a convenient basis for the hy- 
drodynamics of fluids with zero equilibrium energy density, using (p, u") as the fundamental 
variables, and we provide the relations which can be used to express other linearly depen- 
dent fluid scalars (up to third order in gradients), vectors and tensors (up to second order in 
gradients) in terms of this basis. 

At zeroth order, there is only one scalar, p, and one vector, u", along with one symmetric 
tensor orthogonal to n°, hat = rjab + UaUb- At first order in gradients, the scalar daU^, or 
equivalently IC^, is higher order in gradients due to the incompressibility equation. Therefore 
Dlnp is the only independent scalar at this order. The vectors orthogonal to u"" are the ac- 
celeration a" and the pressure gradient orthogonal to the fluid velocity Z?" In p. However, the 
equations of motion imply that only one of them, which we choose to be -D^. is indepen- 
dent. The symmetric tensors orthogonal to u"" are ICab and Dlnphab- The latter is isotropic, 
i.e., proportional to hab- There is an one-to-one mapping between non- isotropic symmetric 
tensors and traceless symmetric tensors. Isotropic symmetric tensors (proportional to hab) 
may be classified by their multiplicative prefactor. Putting together these considerations, 
one can derive a basis for scalars, vectors and tensors up to first order in gradients, which 
is summarised in Table 1. Dependent scalars and vectors can be expressed in terms of this 
basis using the first order equations of motion given in (A.l). Note that derivatives of the 
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Order in gradients 





1 


Scalars 


P 


D Inp 


Vectors orthogonal to u'^ 




Dl Inp 


Symmetric non-isotropic tensors orthogonal to u" 







Table 1. Basis of scalars, vectors and tensors at zeroth and first order in gradients. 



fluid velocity can be expressed as 

daUb = Kab + ^ab " Uattb ■ (A. 2) 

At second order in derivatives, one can use the fluid equations (A.l) and their first 
derivative to derive the following relationships: 

daDl Inp = -JCabJC""' + Qabfi'"'' + 0(8^), (A.3) 

D^Dl Inp = h''''dadb hip + 0{d-^) = -KabK."^ + ^ab^"^ + {D± lnpf + 0{d^), (A.4) 

uVSc^rflnp = DDlnp+{D±lnpf + 0{d^), (A.5) 

hldbi}\ = hldbK\ + {ICab - ^ab)Di Inp + Oid'), (A.6) 

K^^DKcd = DIC = 0{d^), (A.7) 

h^D^^D^ Inp = hlhfdcddlnp + KabDlnp + 0{&^), (A.8) 
d^^aDi) Inp = UaUb{D± \npf + ICabDlnp + hlhfdcddlnp 

+ u^^a{^{Kb)c + nb)c)Dl\np-Di^D\np + Di^\npD\np)+0{^^), (A.9) 

KhiDlCcd = -Khidcdd Inp - JCabD Inp + InpD^ Inp - K^JC^b - K'^cb 

+ 0{d^). (A.IO) 

Note also the following exact relation 

[D, D^] = aaD + Uaa'Di - (/C^ + nJ')D^ . (A.ll) 

Taking these relations into account, we can choose as a basis the five scalars and the five 
vectors orthogonal to the fluid velocity that are indicated in Table 2. One also finds six 
independent symmetric non-isotropic tensor fields orthogonal to u". Dependent quantities 
may be expressed in terms of this basis using the relationships given above. Our result is 
consistent with that of [35]: the number of independent fields at second order coincides (even 
though the basis of scalars, vectors and tensors is different). 

At third order in gradients, we obtain the following relationships using the fluid equations 
(A.l) and their first and second derivatives: 

IC^DICab = -IC^D^Di Inp + IC^D^ InpD^ Inp- K^lCllCl - IC^n^n"^ 

+ 0(a^), (A.12) 
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Order in gradients 


2 


Scalars 


ICablC^', (Dlnpf, DDlnp, (D^lnpf 


Vectors orthogonal to u"" 


Kd^lCl D^Dlnp, ICabD'ilnp, nabD'^lnp, \npD\np 


Symmetric non-isotropic tensors 
orthogonal to u"" 


ICyCcb, /C^^f^idfe), ^a^cb, h^^hfdcddlnp, ICabDlnp, 
In p In p 



Table 2. Basis of scalars, vectors and tensors at second order in gradients. 



Order in gradients 


3 


Scalars 


fC^'D^D^lnp, KlKlKl, lCia\^\, Dlnp^ab^'^', dafC^Dl Inp, 
IC^D^lnpD^lTip, DlnplCabJC'', D^lnp, D'^lnpDlnp, {DlnpY 
{D^ \npfD\ap, DD^hipD^ \np 



Table 3. Basis of scalars at third order in gradients. 



n^'^DQab = -VLab^''^D\Tip + 2K.in\9.\ + 0(5^), (A.13) 
DdaDl Inp = IIC^D^D^ Inp - IIC^D^ InpD^ lnp + 2/C^/C^/C;^ + GlC^n^n^^ 

-2nabn''''Dlnp + 0{d^), (A.14) 

IC^'dadblnp = IC^D^D^lnp- ICablC^Dlnp + 0{d^), (A.15) 

dan^D'i Inp = daK,\D\ Inp + IC'^D^ InpD^ Inp + 0{d^), (A.16) 

n^^'D^D^ Inp = Vlab^^^D Inp + 0(9^), (A.17) 
L»l InpD^D lnp = Dl InpDD^ Inp + {D± In pfD Inp + IC^D^ In pD^ Inp 

+ 0{d^), (A.18) 
daD'lD Inp = AlC^D^D^ Inp - /C'^L'f In pD^ Inp + + QK.I^\Vl\ 
- ICablC^^D Inp- ^ab^^^D Inp + 3L»1 InpDL'f Inp 

+ (L»xlnp)2z)lnp + 2aa/C"feLl^lnp + 0(a^), (A.19) 
D^^dalCl = -A]Cl^\9.\ + dafC^D^ Inp + /C^^Df InpD^ Inp 

+ Sn^b^^^'D Inp - IC'^D^D^ Inp + (9(9^). (A.20) 

Note that the higher order correction terms to the fluid equations (A.l) only appear at 
subleading orders in gradients in these relations. Taking these relations into account, we may 
choose a basis of twelve scalars^ at third order, as indicated in Table 3. 



^Note added: These twelve scalars appear to be consistent with the classification obtained recently in [36]. 
Indeed, according to Tables 4, 8, 9 and 10 of [36], at third order there is one scalar involving three derivatives 
of a zeroth order quantity, four scalars involving the product of a one-derivative term and a two-derivative 
term, and seven scalars representing products of three one-derivative terms. Note however that the results of 
[36] require the use of Landau gauge, and so are not directly applicable to the Rindler fluid with vanishing 
equilibrium energy density considered here. 



- 26 - 



References 



I. Bredberg, C. Keeler, V. Lysov and A. Strominger, "From Navier-Stokes To Einstein", 
[arXiv:1101.2451]. 

G. Compere, P. McFadden, K. Skenderis and M. Taylor, "The Holographic fluid dual to vacuum 
Einstein gravity," JHEP 1107 (2011) 050 [arXiv:1103.3022 [hep-th]]. 

I. Bredberg and A. Strominger, "Black Holes as Incompressible Fluids on the Sphere," 
arXiv: 1106.3084 [hcp-th]. 

T. -Z. Huang, Y. Ling, W. -J. Pan, Y. Tian and X. -N. Wu, "From Petrov-Einstein to 
Navier-Stokes in Spatially Curved Spacetime," JHEP 1110, 079 (2011) [arXiv: 1107. 1464 

[gr-qc]]. 

R. Nakayama, "The Holographic Fluid on the Sphere Dual to the Schwarzschild Black Hole," 
arXiv: 1109. 1185 [hcp-th]. 

D. Anninos, T. Anous, I. Bredberg and G. S. Ng, "Incompressible Fluids of the do Sitter 
Horizon and Beyond," arXiv: 11 10.3792 [hep-th]. 

G. Chirco, C. Eling and S. Liberati, "Higher Curvature Gravity and the Holographic fluid dual 
to flat spacetime," JHEP 1108, 009 (2011) [arXiv: 1105.4482 [hep-th]]. 

C. Niu, Y. Tian, X. -N. Wu and Y. Ling, "Incompressible Navier-Stokes Equation from 
Einstein-Maxwell and Gauss-Bonnet-Maxwell Theories," arXiv: 1107. 1430 [hcp-th]. 

T. Damour, "Quelqucs proprietes mecaniques, electromagnctiqucs, thcrmodynamiques et 
quantiques des trous noirs," These de Doctorat d'Etat, Universite Pierre et Marie Curie, Paris 
VI, 1979. 

T. Damour, "Surface effects in black hole physics," Proceedings of the second Marcel 
Grossmann meeting on General Relativity (1982), Ed. R. Ruffini, North Holland. 

K. S. Thorne, R. H. Price, and D. A. Macdonald, "Black Holes: the Membrane Paradigm," Yale 
University Press, New Haven, USA (1986). 

S. Bhattacharyya, V. E. Hubcny, S. Minwalla and M. Rangamani, "Nonlinear Fluid Dynamics 
from Gravity," JHEP 0802 (2008) 045, [arXiv:0712.2456]. 

I. Fouxon and Y. Oz, "Conformal Field Theory as Microscopic Dynamics of Incompressible 
Eulcr and Navier-Stokes Equations," Phys. Rev. Lett. 101 (2008) 261602, [arXiv:0809.4512]. 

S. Bhattacharyya, S. Minwalla and S. R. Wadia, "The Incompressible Non-Relativistic 
Navier-Stokes Equation from Gravity," JHEP 0908 (2009) 059, [arXiv:0810.1545]. 

C. Eling, I. Fouxon and Y. Oz, "The Incompressible Navier-Stokes Equations From Membrane 
Dynamics," Phys. Lett. B 680 (2009) 496, [arXiv:0905.3638]. 

G. Policastro, D. T. Son and A. O. Starinets, "The Shear viscosity of strongly coupled N=4 
supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87, 081601 (2001) [hcp-th/0104066]. 

R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, "World- Volume Effective Theory for 
Higher-Dimensional Black Holes," Phys. Rev. Lett. 102, 191301 (2009) [arXiv:0902.0427 
[hep-th]]. 

R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, "Blackfolds in Supergravity and String 



-27- 



Theory," JHEP 1108, 154 (2011) [arXiv:1106.4428 [hep-th]]. 

R. -G. Cai, L. Li and Y. -L. Zhang, "Non-Rclativistic Fluid Dual to Asymptotically AdS 
Gravity at Finite Cutoff Surface," JHEP 1107, 027 (2011) [arXiv:1104.3281 [hep-th]]. 

J. Mei, "Towards A Possible Fluid Flow Underlying the Kerr Spacetime," arXiv: 1104.3728 
[hep-th]. 

V. Lysov and A. Strominger, "From Pctrov-Einstein to Navier-Stokes," arXiv: 1104.5502 
[hep-th]. 

S. Kuperstein and A. Mukhopadhyay, "The unconditional RG flow of the relativistic 
holographic fluid," JHEP 1111 (2011) 130 [arXiv:1105.4530]. 

D. K. Brattan, J. Camps, R. Loganayagam and M. Rangamani, "CFT dual of the AdS 
Dirichlet problem : Fluid/Gravity on cut-off surfaces," arXiv:1106.2577 [hep-th]. 

H. Verschelde and V. I. Zakharov, "Notes on relativistic superfluidity and gauge/string 
duality," arXiv:1106.4154 [hep-th]. 

C. Eling and Y. Oz, "Holographic Screens and Transport Coefficients in the Fluid/Gravity 
Correspondence," Phys. Rev. Lett. 107, 201602 (2011) [arXiv:1107.2134 [hep-th]]. 

F. G. Rodrigues, W. A. Rodrigues, Jr and R. da Rocha, "Maxwell and Navier-Stokes Equations 
Equivalent to Einstein Equation," arXiv:1109.5274 [math-ph]. 

T. -Z. Huang, Y. Ling, W. -J. Pan, Y. Tian and X. -N. Wu, "Fluid/Gravity duality with Petrov 
boundary condition in a spacetime with a cosmological constant," arXiv:1111.1576 [hep-th]. 

D. Marolf and M. Rangamani, "Causality and the AdS Dirichlet problem," arXiv:1201.1233 
[hep-th]. 

E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge theories," 
Adv. Theor. Math. Phys. 2 (1998) 505 [arXiv:hep-tli/9803131]. 

S. Bhattacharyya, V. E. Hubeny, R. Loganayagam, G. Mandal, S. Minwalla, T. Morita, 

M. Rangamani and H. S. Reall, "Local Fluid Dynamical Entropy from Gravity," JHEP 0806, 

055 (2008) [arXiv:0803.2526 [hep-th]]. 

C. Ehng, A. Meyer and Y. Oz, "The Relativistic Rindler Hydrodynamics", arXiv:1201.2705. 

I. Papadimitriou and K. Skenderis, "AdS / CFT correspondence and geometry," 
hep-th/0404176. 

I. Booth, M. P. Heller and M. Spalinski, "Black Brane Entropy and Hydrodynamics," Phys. 
Rev. D 83 (2011) 061901 [arXiv:1010.6301 [hep-th]]. 

I. Booth, M. P. Heller, G. Plewa and M. Spalinski, "On the apparent horizon in fluid-gravity 
duality," Phys. Rev. D 83 (2011) 106005 [arXiv:1102.2885 [hep-th]]. 

P. Romatschke, "Relativistic Viscous Fluid Dynamics and Non-Equilibrium Entropy," Class. 
Quant. Grav. 27 (2010) 025006 [arXiv:0906.4787 [hep-th]]. 

S. Bhattacharyya, "Constraints on the second order transport coefhcients of an uncharged 
fluid," arXiv:1201.4654 [hep-th]. 



- 28 - 



